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Abstract 

The conductance- voltage spectrum of molecular nanostructures measured by 
scanning tunneling spectroscopy (STS) is generally assumed to reflect the lo- 
cal density of states of the molecule. This excludes the possibility of observing 
negative differential resistance (NDR). We report here the observation of NDR 
in the scanning tunneling microscope (STM) current-voltage (I-V) character- 
istics of self-assembled monolayer (SAM) of 4-p-Terphenylthiol molecules on 
gold substrate measured using a platinum probe. We argue that the NDR 
arises from narrow structures in the local density of states at the tip apex 
atom and show that depending on the electrostatic potential profile across 

the system, NDR could be observed in one or both bias directions. 
PACS number: 73.61.Ph, 61.16.Ch 
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Electron transport through molecular nanostructures has been widely studied in recent 
years, using either the scanning tunneling microscope or mechanically controllable break 
junctions0@. The conductance- voltage spectrum, dl/dV, is commonly assumed to reflect the 
local density of states (LDOS) of the moleculeS dl/dV ~ p(E = Ep + eV), Ep being the 
equilibrium Fermi energy. A more accurate description of the conductance- voltage spectrum, 
as shown in Ref. [l], is to take a weighted average of the density of states: 

dl/dV ~ i]p(E = /ix) + (1 - v )p{E = /i 2 ) (1) 

Here p,\ and P2 are the electrochemical potentials in the two contacts and the factor 77 
describes the voltage division across the molecule: p,\ = Ep — r/eV and /i 2 = Ep + {l — rj)eV. 
Based on this viewpoint, NDR, i.e., a negative slope in the I-V curve, cannot occur since 
the density of states is nonnegative. Experimentally, however, we often observe negative 
differential resistance for monolayer of long molecules self-assembled on gold substrate such 
as the one shown in Fig. [I]. 

A possible scenario that can lead to NDR at the atomic level was studied theoretically 
by Lang0 and experimentally by Avouris and co-workers@. The essential argument is that 
if there is a weak link between two parts of the conducting system, each of which has 
relatively narrow features in the density of states in the energy range of interest, NDR is 
likely to occur. This is readily understood from the transfer Hamiltonian point of viewi, 
which relates current to the product of the density of states pl and p^ on both sides of the 
weak link. As the bias is changed, current can decrease if two narrow structures in pi and p^ 
move away from alignment. For the structure we studied (Fig. |]), the weak link is the STM 
tip-molecule junction. If the tip has a featureless density of states, as is implicitly assumed 
in the derivation of Eq. ([I]), we wouldn't expect this scenario to apply. However, narrow 
features in the density of states can develop at the tip atoms in a realistic model of tip. 
This has been widely recognized in the STM study of surfaces! and is used to explain the 
NDR in STM I-V characteristics of boron-exposed silicon surfacei. More recently, Yeyati et 
al.0 have studied the electronic structure of a sharp gold tip in the context of conductance 



quantization in gold atomic-size contacts and shown that narrow resonant states can develop 
at the tip apex atom provided the tip geometry is sufficiently sharp. 

In this paper we explain the occurrence of NDR by taking into account the electronic 
structure of sharp platinum tip (used in our STM measurement) as well as the electrostatic 
potential profile across the tip/molecule system. Using a tip model similar to Yeyati et al.@, 
we find narrow structures in the LDOS of tip apex atom below the equilibrium Fermi energy 
Ep (Fig. H(a)). Since the LDOS of molecule also exhibits sharp structures, we expect NDR 
to occur under applied bias. However, since the narrow structures of LDOS at tip apex atom 
are below Ep, NDR will occur only at positive sample bias if the electrostatic potential at 
the tip apex atom is the same as that of the tip support. But NDR can occur in both bias 
directions if a significant amount of the voltage is dropped between the tip apex atom and 
the tip support. We argue that for very sharp tip geometry, this is likely to be true when 
the tip is close to the molecule. Since the tip is composed of a small cluster of platinum 
atoms, the screening length in this region can be much larger than that in the bulk and also 
larger than the cluster size, making it possible to maintain an electrostatic potential drop 
between the tip apex and the tip support under applied bias. Experimentally, we observe 
NDR in both bias directions when the tip is very close to the molecule, in agreement with 
the theory. 

Our model is illustrated in Fig. |2|(a). The self-assembled monolayer (SAM) of 4-p- 
Terphenylthiol molecules are synthesized using standard procedure. I-V data presented here 
(Fig. |l]) represents the average of 25-50 consecutive I(V) sweeps taken at a fixed position on 
the sample!. The SAMs attach strongly to the gold (111) surface through the sulfur end atom 
forming a strong chemical bond with good orbital overlap0. The STM tip usually couples 
weakly to the molecule, corresponding to the physisorption situation. As in our previous 
studies^, we'll use the Extended Hiickel Theory (EHT) to describe the whole molecule-STM 
system (for a recent justification of using EHT in STM study of metal surfaces, see Ref. |TTD , 
taking into account the 5d6s6p orbitals of platinum and gold. 

Similar to Yeyati et al.i, we model the tip geometry as a small cluster of Pt atoms stacked 



on the (111) surface of semi-infinite supportO, originating from a monatomic apex, adding 
3 nearest neighbor atoms to the 2nd layer and 7 nearest neighbor atoms to the 3rd layer 
(we also include 14 nearest neighbor atoms on the surface of the tip support and 6 atoms 
on the gold substrate surface in the calculation of I-V characteristics). The on-site energies 
for the orbitals of the 3 layers of tip atoms are modified self-consistently by adjusting the 
occupation number of each orbital to that of the neutral atom until local charge neutrality 
is achieved on each atomic site0HIHIl. 
We calculate current usingEl: 
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where f{E) is the Fermi distribution and fix, [i2 are the electrochemical potential of the gold 
substrate and tip support respectively, The transmission T(E, V) can be 

calculated using the scattering theory of transport, as we shall describe shortly, but we can 
get more insight if we use the transfer Hamiltonian formalism! to relate the transmission to 
the local density of states on either side of the STM tip-molecule junction: 

T(E, V) = 4vr 2 | M LR | 2 p L (E - eV L )p R (E - eV R ) (3) 

where Mlr is the coupling matrix element and Vl and V R are the electrostatic potentials of 
the molecule and the tip apex atom respectively. 

The LDOS of both the molecule and the tip apex atom show narrow features (Fig. 0(a)), 
which are calculated from the Green's function of molecule-tip system usingi: 

Pmolecule{tip) = ~—Tr{Im(GS) mo i ecu i e ^ip)} (4) 

where G(E) is defined by G(E)(ES-H) = (ES-H)G(E) = /, and S is the overlap matrix. 
The narrow structures of LDOS at tip apex can be understood qualitatively from the fact 
that platinum is third row transition metal, whose LDOS around the equilibrium Fermi 
energy Ep is mostly due to the contribution of the 5d orbitals, located slightly below Ep. 



For the given sharp tip geometry, the coupling of the tip apex atom to its local environment 
is much weaker comparing to that in the bulk, so the levels remain fairly sharp. 

What is the electrostatic potential profile? Applying an external bias changes the relative 
electrochemical potential of the gold substrate and tip support, which are taken as infinite 
electron reservoirs. The electrostatic potential (p(r) is determined self-consistently from the 
Poisson equation V 2 Sif(r) = eSn{r). We only need to calculate the change in electrostatic 
potential 8ip(r) (which in turn modifies the molecular Hamiltonian) since its equilibrium 
value has been included in the equilibrium Hamiltonian. The electrostatic potential change 
in each electrode follows that of the electrochemical potential which provides boundary 
condition to the Poisson equation. The important point is that we take the 3 layers of tip 
atoms and the molecule on an equal footing, viewing the whole molecule-tip system as an 
"extended molecule" sandwiched between the two electrodes. In the region near the tip 
apex, both the electrochemical and the electrostatic potential can be different from those 
inside the tip support and from each otherS. 

As a first approximation, if we neglect any charge buildup within the molecule-tip sys- 
tem, and assume that the gold substrate and tip support act as two infinite parallel plates of 
a capacitor, then the electrostatic potential varies linearly, as shown in Fig. 0(b). For sim- 
plicity, we assume that the molecular energy levels simply float up by an amount around the 
average electrostatic potential change, that is, we neglect any Stark shift of the energy levels 
due to the electric field inside the molecule^. The origin of NDR can be best understood 
by examining how the LDOS structures of the molecule and the tip apex atom sweep past 
each other under applied bias (Fig. ^). NDR will occur for both bias directions if we assume 
a significant voltage drop between the tip apex and the tip support, i.e., eV ape x < eV, as 
evident from Fig. [| 

It can be seen from Fig. |](a) that this approach predicts the I-V characteristics quite 
well with four fitting parameters, namely the tip-molecule distance, the equilibrium Fermi 
energy E F , the constant coupling matrix element M LR and the electrostatic potential change 
of the molecule. We also calculated the transmission function using the equation (see page 



148 in Ref. 0): 

T{E, V) = Tr{T L (E)G R (E, V)T R (E - eV)G A (E, V)} (5) 

which goes beyond the transfer Hamiltonian formalism and takes arbitrary sample/electrode 
coupling into account. Also, a better approximation to the electrostatic potential should 
take charging effects within the molecule into account. In tight-binding theory, this gives 
Stfimoi = Vo + U5n mo i/e, where V is the average elctrostatic potential change of the molecule 
without charge buildup and Sn mo i is the excess number of electrons in the molecule induced 
by the applied bias. Since the molecule is strongly coupled to the gold substrate, we can 
assume that it is in equilibrium with the gold substrate and get Sn mo i by solving the following 
self-consistent equation: 
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The charging energy U can be estimated from the energy gap in molecular optical spectra, 
electron affinity and ionization potential of the molecule: U — (I — A — E g )/2 — Wi m , where 
E„ is the optical gap and Wi m ~ .f is the image potential due to the substrate£3(i? mo ; 
is the distance from the center-of-mass of molecule to the substrate). Since I — A ~ 4eV, 
E g ~ 2.5eV, and R mo i ~ 8.5A, U is ~ 0.3eV. We have used U = O.leV in our calculation. 
The resulting electrostatic potential change is included in the calculation of the Green's 
function and transmission probabilities using Eq. (|5|). Results of such calculation show 
similar qualitative features to that obtained from the simple transfer Hamiltonian argument 
(Fig. 1(b)). 

In conclusion, we have presented a simple explanation of the NDR in STM I-V character- 
istics of self-assembled monolayer of 4-p-Terphenylthiol molecules in terms of the electronic 
structure of sharp platinum tip. The major approximation in our theory is the use of the 
Extended Hiickel theory and a simplified treatment of electrostatic potential variation. A 
central factor in understanding the NDR is the narrow structures in the LDOS at the tip 
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apex which has been noted by other authors in different contexts. However, this alone is not 
enough to explain the occurrence of NDR in both bias directions. NDR for negative sample 
bias can be understood only if we allow for an electrostatic potential drop between the tip 
apex and the tip support. This is possible for very sharp tip geometries, where screening 
length is large. But if the tip is relatively flat, the potential drop between the tip apex atom 
and the tip support will be reduced. Indeed, experimentally we do observe NDR only in 
positive sample bias direction for some of our tips. However, even with these tips, we ob- 
serve NDR in both bias directions when we move the tips closer to the molecule, indicating 
a significant increase in the potential drop between the tip apex and the tip support. It's 
not clear to us why the potential drop should increase significantly, since the distance range 
by which the tips move is relatively narrow. Perhaps a more complete theory is needed that 
involves self-consistent evaluation of the electronic structure for the entire electrode/sample 
system, such as described by LangEl. Our main purpose here, is then to show that the elec- 
trostatic potential profiles play a crucial role in determining the shape of I-V characteristics 
of molecular nanostructures — a role that has not been adequately recognized. 

This work was partially funded by the National Science Foundation under Grant No. 
9708107-DMR. One of us (R.R.) would like to thank Norton Lang for a stimulating discussion 
during the very early stages of this work. 
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FIGURES 




4-p-Terphenylthiol 
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Sample Voltage (V) 

FIG. 1. Experimental I-V characteristics of the molecule as a function of the molecule-tip 

distance. A negative differential resistance is observed for small tip/SAM separation. The set 
point voltage is -5V; the set point currents range between 3.05 nA and 30.5 nA. 
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(b) 



Tip Support 




Gold Substrate 

FIG. 2. (a) SAM of 4-p-Terphenylthiol on gold substrate (only one molecule is shown here for 
clarity), also shown is the STM tip. (b) Electrostatic potential profile under applied bias, here we 
only show the case of positive sample bias. 
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At Equilibrium 
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FIG. 3. Illustration of how the LDOS structures of molecule and tip apex atom sweep past each 
other under applied bias. The electrochemical potential of gold substrate \x\ is taken as energy 
reference (The LDOS curves have been horizontally offset for clarity). Parameters used are the 
same as those used for the I-V calculation in Fig. || (a), (a) At equilibrium, \i\ = //2 = Ep; (b) 
At negative sample bias, the LDOS curves of molecule and tip apex atom float down relative to 
gold substrate; (c) At positive sample bias, the LDOS curves of molecule and tip apex atom float 
up relative to gold substrate. NDR already occurs before sample bias reaches ±3.5(V). 
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FIG. 4. (a) I-V characteristics calculated using transfer Hamiltonian theory for tip-molecule 
distance= 3.50A. Parameters used are: [ M LR [ 2 = 8.2 x l(T 8 (ey) 2 , E F = -11.05(eV) and the 
electrostatic potential change of the molecule is taken as 1.15 times the average electrostatic po- 
tential change, (b) I-V characteristics calculated using Eq. (^) for various tip-molecule distance d 
(in unit of angstrom), here we use Ep = —11.15 (eV), U = O.l(eV). 
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